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Abstract 

We investigate the thermodynamic properties of the nonlocal tachyon in 
string field theory. We use previously developed perturbative methods for non- 
local fields to calculate the partition function and the equation of state in the 
high temperature limit. We find that the tachyon condensation can be seen as 
a second order phase transition. We compare our results with those of ordinary 
scalar field theory. We also calculate the one loop finite temperature effective 
potential. 
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1 Introduction 



Tachyons are ubiquitous in string theory. They made their first appearance in closed 
bosonic string theory. Since then, they have been found to exist as open string 
excitations in the world volume of D-branes in bosonic theories and in non-BPS D- 
branes in superstring theory, as well as excitations in brane-anti-brane systems; see pQ 
for a review. In conventional field theory the appearance of a tachyon usually implies 
that we are perturbing around an unstable vacuum; the tachyon should evolve to its 
true vacuum with positive mass-squared around which we can perform perturbative 
quantum calculations meaningfully. There is growing evidence to suggest the same is 
true in string theory [TJ |2] . For instance, the open string tachyons are simply thought 
to represent the instability in the various D-brane systems. The rolling of the tachyon 
from the unstable potential hill to its stable minimum then describes the dissolution of 
the unstable D-branes into closed string excitations around the true vacuum which no 
longer supports open string excitations. This process is often referred to as tachyon 
condensation in the string literature. While the classical dynamics of the above 
process has been studied extensively, relatively little attention have been paid to the 
quantum theory of the tachyon. 

The present paper will demonstrate that, just as in ordinary field theory, one can 
consistently perform quantum calculations by perturbing around the true tachyon 
vacuum in string theory. For this purpose, we will consider the simplest level of 
truncation in the String Field Theory (SFT) approach |3] where only the tachyon field 
is kept [I]. Work done on such string-inspired nonlocal models have been fruitful in 
the context of string theory [5]- [TO], cosmology [U]-[I7j and particle physics [IS] . (For 
related work on nonlocal gravity see [E]-[26].) With regard to quantum phenomenon 
in string theory these simpler models, for instance p-adic strings [27] , have been shown 
to reproduce properties such as thermal-duality [2H] (which has been variously argued 
in the string literature [2H]) and Regge behavior [30]-[3T]. 

More specifically, in this paper we will study the possibility that tachyon con- 
densation can be seen as a phase transition at high temperature. Schematically, the 
string field \l/ can be thought of as a matrix- valued 1-form [I] with a Chern- Simons 
type action. For instance, the bosonic SFT action is given by 



where Q is the BRST operator which is normalized to provide canonical kinetic terms 
to the various particle fields contained in the SFT spectrum, g is the string coupling 
constant, and l/y/a' is the string tension. The * product has the effect of diffusing 
or smearing out the interaction over the string length. For instance, if / corresponds 
to a canonically normalized particle field in the SFT spectrum, then an interaction 
term involving the * product in the string field \1/ translates into interactions for / 
where they only enter in momentum-dependent combinations 





/ = exp a 
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Equivalently, one can work with the redefined fields / in terms of which the interac- 
tions have the usual polynomial form, but the kinetic operator picks up the nonlocal 
exponential derivative dependence [I]. Thus, if we keep only the tachyon field, then 
the corresponding field theory action has the form 



where m 2 is the mass of the tachyon at the maximum, and M is the scale of nonlocality 
that describes stringy interactions. Both m and M are proportional to the string 
tension. For example, for the bosonic SFT m 2 = — 1/a', while M 2 ~ 1/a! ln(3\/3/4) . 
We will treat m and M as independent parameters for technical and phenomenological 
reasons. The V(<ft) represents a polynomial interaction, typically cubic or quartic. 

In this paper we employ finite temperature methods that were developed for 
such nonlocal theories in [32j [28] to investigate the thermodynamic properties of 
the tachyonic excitations. We will discover that, just as in ordinary Quantum Field 
Theory (QFT), the SFT tachyon undergoes a second order phase transition. The 
effect of the stringy nonlocality seems to weaken the phase transition. This means 
that the discontinuity in the specific heat as one approaches the critical temperature 
from above and below decreases as M decreases. We emphasize that the QFT limit 
is expected to be attained in the limit M — > oo, and we explicitly verify that this is 
indeed the case. In this paper we work in the limit M 3> m; that is, we are close to 
the particle limit, but the formalism and techniques that we have developed can be 
employed to understand the (perhaps) more interesting and relevant situation where 
M ~ m. This will be reported in a future paper. 

Our calculations help to clarify the relation between the conventional QFT, which 
is based on a renormalization prescription, and the string inspired nonlocal actions 
where loop diagrams are typically finite. Intuitively, the exponential cut-off scale, 
M 2 , acts as a Lorentz-invariant physical regulator. The expressions for the various 
thermodynamic quantities in these nonlocal models are almost identical to what one 
obtains using traditional renormalization prescription, except that there are correc- 
tions which are suppressed as 0[exp(— M 2 /4T 2 )]. This happens because one takes 
the limit M — > oo after imposing the "renormalization conditions" , according to the 
standard renormalization prescription, whereas the M in stringy Lagrangians is a 
finite physical parameter encoding the nonlocality of the model. This opens up pos- 
sibilities for phenomenological applications in particle physics if M is close enough to 
the scale of Standard Model physics. For one proposed alternative to the Standard 
Model using SFT see [E]. 

Another goal of this paper is to pave the way to explore possible connections be- 
tween string theory and cosmology. In recent years, string thermodynamics has found 
several applications in the early Universe cosmology [33J. In particular, there have 
been efforts to see whether stringy thermal fluctuations can play a role in the forma- 
tion of the anisotropics in the Cosmic Microwave Background Radiation (CMBR). 
There are also the so-called warm inflationary scenarios where particles are continu- 




(1.2) 
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ously produced and which then thermalize and influence inflationary dynamics, both 
at the level of the background and the fluctuations [M] . It would be interesting to 
consider similar scenarios where stringy excitations are produced instead, potentially 
providing a prospect to observe stringy properties in CMBR. Our calculations would 
be relevant for such a discussion. 

The paper is organized as follows: In section 2, we introduce the SFT model 
and its finite temperature formulation. By calculating the temperature dependence 
of the effective mass and minimum of the potential, we demonstrate that a second 
order phase transition occurs in SFT at high temperature, and estimate the critical 
temperature. In section 3, we follow the traditional perturbation theory approach 
to calculate the 1 and 2-loop diagrams contributing to the partition function of the 
tachyon. This enables us to obtain the equation of state for the thermal excitations 
of the tachyons around the true minimum. The same techniques can be used to 
calculate N-point diagrams for arbitrary values of N at high temperature. This is 
sufficient for us to obtain the critical temperature and determine the nature of the 
phase transition. In section 4, we summarize the analytical computations and provide 
numerical results. This enables us to compare the SFT equation of state with the 
analogous QFT equation of state. In section 5, we calculate the 1-loop effective 
potential which extrapolates away from the equilibrium states. Concluding remarks 
are contained in section 6. 



2 Action and Critical Temperature 

In this section we introduce the SFT action and show that a second order phase 
transition is to be expected at high temperature. More elaborate calculations of the 
equation of state and of the effective potential follow in sections 3, 4 and 5. 

2.1 Action 

Our starting point is the SFT type action given by [Hj 



The metric is such that □ = — d 2 + V 2 . In the bosonic cubic string field theory the 
action for the tachyon is of the above form with V(<f)) oc <ft 3 , while for the supersym- 
metric case one expects a quartic potential. In either case the mass-squared term is 
always negative, m 2 < 0, indicating the presence of a tachyon at = 0. In this paper 
we focus on the quartic coupling 



In SFT A is proportional to g 2 . Since the potential is bounded from below, we expect 
to be able to perform loop calculations without encountering any pathologies. In this 





V{ct>) = A0 4 
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context we note that although the presence of higher derivative terms usually indicate 
the existence of additional ghost-like states, the fact that in the SFT action they 
combine into an exponential ensures that there are no extra poles in the propagator. 
Thus there are, in fact, no additional physical states, ghosts or otherwise. There are 
also strong arguments to suggest that the initial value problem in these models are 
well defined, and classical trajectories can be uniquely specified by only a finite set 
of parameters 



There are two interesting limits of the action expressed by eq. (2.1). When m is 
fixed and M — > oo, one recovers the conventional local field theory action for a scalar 
field. When M is fixed and m — > oo, one recovers the p-adic field theory [27 \. In this 
paper we consider the case when both m and M are finite but with \m\ -C M. 

Application of the usual finite temperature formalism to nonlocal actions involv- 



ing an infinite series of higher derivative terms, such as (2.1), have been studied 
recently [28J. The basic prescription is rather straightforward and resembles the fi- 
nite temperature methods implemented in usual local quantum field theories. The 
main difference is that the propagator gets modified by the presence of the nonlo- 
cal terms, and ultraviolet divergences in the quantum loops are either softened or 
eliminated altogether. In some sense the mass parameter M acts as an ultraviolet 
cutoff. Thus the renormalization prescription is somewhat modified from the usual 
field theories. 

One of the important differences between SFT and p-adic theory is that, in the 
latter case, = corresponds to a minimum while in the SFT case it is a maximum, 
hence perturbative calculations around = are not well-defined. (This is elaborated 
on in more detail in the following sections.) To perform quantum loop calculations we 
must expand around the true minimum. At the classical, or tree, level the minimum is 
located at O = /x/2\/A, where we have defined /x 2 = — m 2 > 0. See Fig. 1. However, 
at finite temperature the minimum is shifted to smaller values of 0. To account for 
this we expand around the true minimum, v(T): 

= v(T) + <P f (2.3) 

where v(T) is independent of space and time but does depend on the temperature, 
while 0/ is the fluctuation around it whose average value is zero. When doing loop 
calculations using the fluctuation 0^ around the average value of the field v(T) at 
the temperature T, we should use the dressed propagator (or our best estimate of 
it) and not the bare propagator. As in conventional models of spontaneously broken 
symmetries in local field theory, we will see that as the temperature increases the 
values of the condensate v(T) and the effective mass both decrease to zero at a critical 
temperature T c . This can be done in the usual way by introducing an additional 
mass-shift 5m 2 in the the quadratic part of the action and then subtracting it as a 
counter-term. Then the Lagrangian reads 

£ = £ quad + Ant + Cct - V c i(v) (2.4) 
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where 




Figure 1: The vacuum potential V in dimensionless form. 



•^quad 

V cl (v) 



\<f> f [e- u/M2 (□ + /i 2 ) - I2\v 2 - Sm 2 \ fa 
\bm 2 ^) 



(2.5) 
(2.6) 



The v(T) and Sm 2 (T) are as yet undetermined functions of the temperature. 



2.2 Critical Temperature 

The functions v(T) and 5m 2 (T) introduced above can be determined by the 1-point 
and 2-point loop diagrams, respectively. In this section we derive formulas for them 
at the 1-loop mean-field level. We begin by noting that the thermodynamic potential 
is given by 



n = vjv 



T 
V 



In 



[<%] exp ( / dr I d 3 x [£ quad 
o Jv 



ct 



(2.7) 



where /3 — 1/T and V is the volume. The thermodynamic potential must be a 
minimum with respect to variations in v, thus providing an equation to determine 
v(T). Similarly the function Sm 2 is determined by the Schwinger-Dyson equation in 
the 1-loop mean-field approximation. 

According to the formalism that was developed in [28J to compute Feynman 



diagrams in nonlocal field theories such as (2.1), the only diagrammatic rule that 
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needs to change is the propagator. From (2.5) we see that 



V e ( w a+pa)/Afa ( w 2 + p 2 _ ^2) + 12 A W 2 + £ m 2 ( 2 ' 8 ) 

where we have introduced the Matsubara frequencies, u n = 2nirT, as appropriate for 
finite temperature field theory [36J. We emphasize the exponential suppression of the 
propagator at large momenta which is a typical characteristic of these theories and 
which helps to regulate the ultarviolet divergences of loop diagrams. 

Let us now try to compute these quantities at the 1-loop level; a more rigorous 
comprehensive analysis will be provided in the next two sections. 

*| = => (i\v*-S)v + 12\vT W** C = (2.9) 

n ^ 

where the first term comes from the classical potential while the second term comes 
from the 1-loop tadpole diagram. For sufficiently small T there is a local maximum 
at v = and a minimum at 



3T E / ^V(u n , P] v,5m 2 ) (2.10) 



2 = fj?_ orT1 ^ f d 3 p 
' " 4A 



At a critical temperature T c these become degenerate. The value of T c is determined 
by 

^2 f - d3p 



T = 
c 12A 



(2tt) : 



(2.11) 



so that above this temperature there is only one minimum at v — 0. This temperature 
will later be identified with the critical temperature of a second order phase transition. 

The two equations above are not sufficient to determine v(T) and T c because the 
right side depends on T> which itself depends on 5m 2 . In the mean-field approximation 
it is determined by the 1-loop self-energy diagram which involves only the 4-point 
vertex, not the other diagram involving the 3-point vertices. (For details see the next 
section.) This is a gap equation. 

d 3 p 



5m 2 = 12XTJ2 J -^V(u ni p;v,5m 2 ) (2.12) 



Equations ( 2.10[ ) and ( |2.12[ ) are to be solved simultaneously and self-consistently to 
determine v (T) and 5m 2 (T). 

First consider T > T c : the minimum is located at v — and so 5m 2 can be 
determined as the solution to the single equation 



5m z 



12XT J2j ( 2v r)3 e K+ P W (w 2 +p 2_^ + 6m 2 ( 2 - 13 ) 
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For T < T c , on the other hand, there is a simple relation between 5m 2 and v, 
namely 

Sm 2 = ^ 2 - 4Aw 2 (2.14) 

The temperature T c where v goes to zero is exactly the same temperature where the 
effective mass 

m 2 cS = 5m 2 -fx 2 + I2\v 2 = 8\v 2 (2.15) 
goes to zero. The condensate can be determined via 

V = 4A _ ^ J (2tt) 3 e (^+P 2 )/ M2 {u 2 + p 2 - fi 2 ) + fi 2 + 8Xv 2 ^ ' ' 

and this in turn allows for the direct algebraic determination of Sm 2 . 

To obtain simple analytic results, for the moment we shall assume that /i <C M 
and that M <C T c . These assumptions can of course be relaxed albeit at the expense of 
numerical calculations, and they will be presented in the following sections. Therefore 
we focus on temperatures T ^> M. In this situation a nonzero Matsubara frequency 
will contribute an amount which is suppressed by a factor exp(— 47r 2 T 2 /M 2 ) and is 
totally ignorable. Thus 

T> /— --T?^ -=. 2.17) 

^ J (2tt) 3 Att^/tt v ; 

Hence both the condensate v and the effective mass m 2 ff = 12 Xv 2 — fi 2 + 5m 2 must 
vanish at the same critical temperature given by 

VTv^F fi 2 

Tc ~^XM- (2 ' 18) 

This is strongly indicative of a second order phase transition. For consistency we 
need T c 3> M. This results in the limit 

AM 2 «/i 2 «M 2 . (2.19) 



To summarize: in this limit 



v 2 = 



5m 



fj^ _ 3MT 

4A 47rv / 7r 
3XMT 



below T r and 



2 2 QXMT , 2 / 
m 2 eS = 2/i 2 p- = 8Av 2 (2.20) 

7T V 7T 



v 2 = 

r 2 3AMT 
om = 



2 3XMT 2 

Try 7T 



above T c . 
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3 Equation of State 



In this section we perform a more sophisticated analysis of the equation of state. 
We proceed analytically as far as possible and defer numerical calculations to a later 
section. Readers primarily interested in the results may skip this section. 



3.1 Formalism 

The Lagrangian we will work with is 

£ quad = \<t> e~ U/M2 (□ + /i 2 ) + W - A0 4 (3.1) 

Here we have added a counter-term ^70 2 with a coefficient 7 which will be adjusted so 
that the value of the condensate in the vacuum is the same as the classical expression 
/x 2 /4A. It also insures that no new poles are introduced into the propagator. As 
before we represent the field in the form 

<P = v{T) + <P f (3.2) 

where v(T) is the equilibrium value of the condensate at temperature T and 0/ is 
the fluctuation around it whose average value is zero. After making this shift, and 
acknowledging that terms linear in <pf will average to zero in the functional integral, 
the Lagrangian can be written as 



where 



£ = £ quad + Ant + Cct ~ Vd(v) (3.3) 

C quad = \<t> } [e- u ' M2 (□ + - I2\v 2 - 5m 2 ] <f> f 

£ int = -4Au$-A# 

Cct = |(<Jm 2 + 7 )^ 

V c i{v) = -i-(v 2 + 1 )v 2 + \v 4 . (3.4) 

An additional mass shift 5m 2 has been added to the quadratic part of the action and 
then subtracted as a counter-term. The reason is that as the temperature increases 
the value of the condensate v decreases to zero at a critical T c , and therefore at some 
temperature below T c the effective squared mass in the propagator becomes negative. 
This just means that we should do our calculations with the dressed propagator (or 
our best estimate of it) and not the mean field propagator. The value of 5m 2 has to 
be determined at each temperature self-consistently just like v does. The V c \(v) is the 
contribution to the classical potential from the condensate field. 

The Feynman rules corresponding to the above action are as follows. The dressed 
propagator V in the imaginary time formalism is given by 
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while the mean field propagator is given by 

p-l = e (rt+P 2 )/M* ^2 +p 2_ ^ + UXv 2 (3 g) 

The thermodynamic potential is 

A = V CI (V) - | In |y [<%] exp Qf rfr ^ <* 3 X t^quad + Ant + £ ct A } (3.7) 

where = 1/T and V is the volume. In a diagrammatic expansion the field 0/ is 
represented by a solid line. The vertices can easily be read off from the expressions 
above. The quartic interaction 0j has the vertex —A and the cubic interaction 0^ has 
the vertex —4Xv. A cross or X represents the counter-term —(5m 2 + 7). 

The thermodynamic potential can be considered a function of the equilibrium 
condensate and a functional of the dressed propagator [37J 135] . 

n = v cl (v)-lTj2 [ ^{^(T^-vv^ + ij + f^n^v) (3.8) 

n J 1=2 

The Qi is the /-loop contribution to the potential. (A counter-term counts as one 
loop in this context.) Extremizing with respect to v removes tadpole diagrams, and 
extremizing with respect to T> removes one particle irreducible diagrams. In the 
same way one could remove two particle irreducible diagrams by introducing dressed 
vertices. The equations that determine the equilibrium solution are 



_ = + l)v + 4Au + 12XVT J2]^V + j; ->±-> = (3.9) 



and 

V-'-V-' = 2j2 6ni( 5 V v V) (3.10) 
1=2 

the latter being the Schwinger-Dyson equation. Due to the functional derivative 
the difference of the inverse propagators is frequency and momentum dependent in 
general. Thus 5m 2 should in principle be the self-energy U(u n ,p). However, in 
the approximations used in this paper a constant 5m 2 will suffice. Terminating the 
expansion at two loops results in the diagrams displayed in Fig. 2. We would like to 
point out that all loop diagrams are UV finite on account of the exponential damping 
of the propagator, and they are IR finite except at a second order critical temperature 
where the mass vanishes (correlation length diverges). 

The behavior we expect is that the zero temperature condensate decreases with 
increasing temperature due to thermal fluctuations. If it goes to zero at a finite 
temperature T c then a phase transition ought to have occurred. Now suppose that 
perturbation theory can be applied, at least if we are not too close to T c . To 1-loop 
order the equation for v(T) is 

2 A* 2 + 7 w ^ [ 
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Figure 2: Two-loop contribution to In Z including combinatoric factors and the 
counter-term. 




Figure 3: One loop contribution to the self-energy including combinatoric factors and 
the counter-term. 



The mass shift Sm 2 can also be calculated. At 1-loop order it receives contributions 
from the diagrams shown in Fig. 3, which are obtained from those in Fig. 2. We will 
neglect the diagram involving the three point vertex; this will be justified a posteriori. 

W-mrW^p--, (3.12) 

Therefore, to 1-loop order 

Sm 2 = /i 2 - 4Xv 2 , (3.13) 
and the propagator to this order is 

V -l = e K+P 2 )/A^ ( w 2 + p 2 _ ^ +fl 2 + gXv 2 (314) 

The temperature T c where v goes to zero is the same temperature where the mass 
5m 2 — /i 2 goes to zero. 

The counter-term coefficient 7 can readily be determined by calculating 5m 2 at 
T = 0. In the T ^ limit 



n 

Under the assumption that fi M the integral is trivial, and the requirement that 
5m 2 = requires 

3AM 2 

7 = -^r (3-15) 

Working to higher order in the loop expansion would give the expansion of 7 in a 
power series in A. 



2tt 
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Now let us consider the contribution of the 2-loop contributions to the thermo- 
dynamic potential. Referring to Fig. 1 they are 



^quartic 3 A 

and 



2 



(3.16) 



ficubk = -48AVT J J^ T Y1 J ^k V ( UJ n 1 ,Pi)'D{u n2 ,p2)V(u ni+n2 ,p 1 + p2) 

(3.17) 

in an obvious notation. As emphasized before, the zero Matsubara mode dominates 
at high temperature, T 3> M, but still T < T c so that the cubic interaction does not 
vanish. (If it did then clearly f2 cu bi C = 0.) Hence 

^quartic ~ AM 2 T 2 (3.18) 

and 

^cubic ~ AVT 2 In (M 2 /\v 2 ) . (3.19) 

The diagram with the cubic vertices is smaller by a factor of Xv 2 /M 2 which tends to 
zero at T c even when including the logarithmic factor. 

Similar conclusions can be reached for the self-energy. Using the dressed propa- 
gator of eq. (2.8) as reference we find that 

2 6AMT 

n quartic (u; n ,p) = 3/i pr- (3.20) 

7T\/7r 



which is frequency and momentum independent and has the limit 2/x 2 at T c . In 
contrast, at large frequency and momentum 

n cubic (w n , p) ~ A I MT 2 exp [- (u 2 n + p 2 ) /M 2 ] (3.21) 

OJn + P Z L V J 

which is exponentially suppressed. At zero frequency and momentum 

„ , . 144A 2 T 2 f°° dkk 2 36AVT 9^2, s/2 m 

U cM = -— I w ^ = --— = -—X^. (3.22) 

which goes to zero at T c . This justifies the neglect of the diagrams involving the cubic 
interactions at high temperature. 

Let us summarize the approximate expression for the thermodynamic potential. 
It includes the two loop diagram involving the four point vertex but not the two loop 
diagram involving the three point vertex; as argued previously, the latter is suppressed 
when A 1 and fi <C M. 

a = - J p + 7 ) v> + \v< + FY. f 7S5 ln OS"*" 1 ) 
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(2^ 



The condensate satifies 



2 /^ 2 + 7 



<i 3 p 



below T c while t> = above T c . The propagator is 



Ut; 2 



below T c and 
above T c , where 



3 (^ +p 2)/M 2 



(w 2 +p 2 -/i 2 ) +5m 2 
d 3 p 



n v 7 

It is useful to express the propagator as 

XT 1 ( Wn , p; m 2 ff ) = e ^>/^ (c 2 + p 2 - /, 2 ) + /, 2 + m 



(3.23) 
(3.24) 

(3.25) 
(3.26) 
(3.27) 

(3.28) 



The quantity m c g is approximately the pole mass when it is small compared to M. 
Rigorously speaking it is the screening mass. At zero temperature m 2 ff = 2/z 2 and 
is always non-negative, vanishing only at T c . Below T c we must solve for v(T) self- 
consistently, and the result then also determines m e g(T). Above T c we must solve for 
Sm 2 (T) self-consistently, and this determines m e g(T). 



3.2 Sums and Integrals 

Let us calculate the relevant sums and integrals. We are interested in temperatures 
T > n (previously, in section 2, we considered T > M) and mass scales M > \x. Let 
us start with the oft-appearing quantity 

(2vr)3 

Note that it is convergent in both the IR and UV. To the desired order 

^Ej^^EJ^^r- (3.29) 

The trick is to use the integral representation 

-(ul+p 2 )/M 2 i poo 

E^^^E/ **^«**»" (3-30) 

n n 
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and then use the function 



oo 



which appears so often in the study of the p-adic theory at finite temperature [28] . 
After integrating over momentum 

/ y e = J}^Lf(r/M) (3.32) 

7 (2vr) 3 w2 +p 2 47rv^F M 7 7 V ; 

where we have defined the function 



Therefore, to the desired accuracy 

(2tt) 3 ~ Airy/n 



T ^f^=^=f(T/M). (3.34) 



The integral over a can be performed numerically, but it can also be calculated 
in the low and high temperature limits. These calculations are facilitated by an 
interesting property of the ( function that 




(3.35) 



When s > y/n 

C(s) = 1 + 2e" s2 + 2e" 4s2 + • ■ • (3.36) 

and when s < y/n 

£(s) = ^ (l + 2e~ w2/s2 + 2e- 4 " 2/s2 + •■•). (3.37) 

Using the above, we can calculate the behavior for T < T and for T > T where 
T = M/2y/7T. For T > T 

/(T/M) = l + 2 f:{e-^ 2 ^ 2 - 2 ^ 

n=l ^ 

where $ is the probability integral 

= -= / dte~ t2 . (3.39) 
V 71 " Jo 

This has the high temperature expansion 

f(T/M) = 1 + J^e"- 2 - 2 /- 2 + O (^e-- 2 - 2 /- 2 ) . (3.40) 
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/2t™T\ 



(3.38) 



For T < T the integral can be broken up into two pieces, one from a = 1 to 
a = «o and another from a = ao to a = oo, where «o = M/2y/7rT. 



S (TIM)^ + » 



M 



4^T 



1 



4ttT 5 



+ 



4y^r 
M 



n=l 



2t™ 2 



The low temperature limit is 



f(T/M) 



— e 



M 



-n 2 M 2 /4T 2 



+ e"™ 2 " - nn (l - $(71^)) 



7T\/7T T 



3 M 







T 



-M 2 /4T 2 



where we have used 



1 + 2 E [ + 2 ) e-nV - 2 ™ C 1 - $ (^)) 

n=l L V / 



7T 



1.07163... 



(3.41] 



(3.42) 



(3.43) 



We have not proven the equality stated above, nor have we found reference to it in 
the literature, but it is true to any numerical accuracy that we have done. 

Finally let us turn our attention to the one loop contribution 

d 3 p 



(2tt)' 



In (^V- 1 ) 



It can be expressed as 



da 2 



d 3 p 

Wf 1 Jo e(^+P 2 )/ M2 ( W 2 + p 2 _ ^2) + ^2 + a 2 



d 3 p 
(2^p 



In 



^2 / e (^+p 2 )/A/ 



In the limit of small /i in comparison to M and T the second term can be written as 



+ ln [P (U n + P )J 



1 _ e -K+p2)/M2) 



(27T) 



M 2 



The first term in curly brackets appears in the p-adic limit and it is zero (28]. The 
second term contributes one massless bosonic degree of freedom. 



d 3 p 



(2tt) 3 

The third term can be written as 



\n[p 2 (oo 2 n +p< 



71 

90 



T + vacuum 



(3.44) 



d 3 p 1 1 



(27r) 3 pe^ - 1 



d 4 p 1 /j, 2 MT 



(2vr)V 



57T\/7r 



f(T/M) 
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+ - — 1 =j(T/M) + vacuum. 
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See [36J for the integrals. Hence, to the desired order 
d 3 p - 2 



(3.45) 



I^E / (^3 111 (/? 2 ^ X K,p; meff)) = -^T 4 - + V ^^=-MTf{T/M) . 

(3.46) 



3.3 Equation of State for M, T > /i 



Now we assemble what we have learned in the limit that M ^> fi and T ^> fi. The 
equation of state is expressed as pressure P{T) = — fl(T) as a function of temperature 
T. The pressure is normalized to zero at zero temperature. 

For T < T c the effective mass and condensate are given as functions of tempera- 
ture by 

3AM 2 



m 



eft 



iXv 2 = 2fi 2 



2tt 2 



(3.47) 



The pressure is 



7T~ 



/i 2 T 2 3/x 2 M 2 



P = — T 4 + 

90 24 



32tt 2 



T f T 



3AM 4 

1287T 4 



T / T 
4 ^M ; M 



n 2 



The entropy density s(T) = dP(T)/dT can easily be computed by using the formula 

MT/M) 



T- 



dT 



f{T/M) - C(2vrT/M) 



(3.48) 
)rmula 

(3.49) 



It is 



. . 2vr 2 /i 2 T 

s(T) = T + - — 

v ; 45 12 



2/ 



T 

M 



3AM 3 



16tt 3 



'a ^ T J T 
40T f — 



3/i 2 M 



57TV7T 



(3.50) 



The energy density is e(T) = —P(T) + Ts(T). 



For T > T c the condensate is zero and the effective mass is determined by the 
formula 

_ V M V M / 



771 



(•ir 



4^ 



(3.51) 



The pressure is 



„ 7T 4 /i 2 T 2 3AM 4 
P = — T + c - 

90 24 256vr 4 



~M V M 



-i 2 



16A 



(3.52) 
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and the entropy density is 



2tt 5 



3AM 3 



s(T) - ^ T 3 + ^ 327r3 ^_ 



V 



T 

M 



/2ttT 



M 



T / T 

4 V s 



- 1 
(3.53) 



Both the condensate and the effective mass vanish at the critical temperature T c 
determined by 



f 



( 



M \M) 4^ 
At this temperature both the pressure 



1 -K^/nji 2 
' 3AM 2 



7T 



u 2 T 2 

P(T) = — T 4 + — - 
[ c) 90 c + 24 



12A 



and the entropy density 



frT1 , 2vr 2 „ /i 2 T c /i 2 M 



( 2nT c 
V M 



2/ 



T 

A c 

M 



(3.54) 



(3.55) 



(3.56) 



are continuous, but the heat capacity cy(T) = Tds{T) / dT is not. Hence this is a 
second order phase transition. 



4 Numerical Results and Comparison to Local 
Field Theory 

The string field theory under study has three parameters: M, [i and A. What matters 
for the equation of state is not absolute magnitudes but relative magnitudes. The 
perturbative analysis we have used requires that A 1. However, since fi 2 > with 
the consequence of spontaneous symmetry breaking in the vacuum, the limit A = 
is not allowed since the theory would not be well-defined. We have assumed that the 
string scale M is large in comparison to the mass scale /i. We have also assumed that 
T is large compared to fi but made no assumption about the ratio T/M. From the 
point of view of conventional local field theory, M acts as an ultraviolet regulator. 
If one takes M — > oo in the action, one recovers the normal 4 field theory with 
spontaneous symmetry breaking. Let us examine the limits T c ^ M and T c ^> M 
analytically before turning to numerical calculations. 

Suppose that M — > oo with /i 2 /A held fixed. Then it is easy to show that 

2 f 2^ 2 (1 -T 2 /T 2 ) if T<T C 

m eS {i ) - | ^ 2 (T2/T 2 _ x) \{T >T C l4 ' ij 

with T 2 = ji 2 /\. The condensate is determined by 8\v 2 (T) = m 2 s {T) when T <T C 
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while v(T) = when T > T c . The equation of state below T c is 

P{T) = I 1 ) T — 



90 24 y 12 

V^ 2 A\ , /i 2 T 
S (T) = 4 I — + — I T — 



90 24 y 6 

,'tt 2 A\ 4 /i 2 T 2 
e (T) = 3 ( — + — ) T — 



and above T c is 



90 24 J 12 

c y (T) = 12(- + -\t*-^- (4.2) 
KV ; V 90 24/ 6 V ' 



P(T) = ,^_AW" 2r2 " 4 



90 48/ 24 16A 

Vvr 2 A \ . /i 2 T 
S ( T ) = 4 ™ - ^ T 



90 48 / 12 

,'tt 2 A\ 4 /x 2 T 2 /i 4 
e ( T ) = 3 ™ - 7^ T + oT~ + 



90 48 / 24 16A 

Cv (T) = i 2 (---)t 3 + ^- (4.3) 

(Corrections to these formulas for large but finite M are suppressed by the factor 
exp(— M 2 /4T 2 ).) Clearly P, s, and e are continuous at T c but cy is not. These 
are well-known, conventional finite temperature field theory results [36J. Because A 
is required to be small, and T ^> /i, the equation of state to first approximation is 
e = 3P. To focus on the effect of interactions, especially near T c , it is useful to define 
the dimensionless interation measure (e — 3P)//i 2 T 2 . From the above 

' if T < T r 



/i 2 T 2 



3P | a 11 J - ic 

if T > T, 



\ T 2 ! (4.4) 



4T 2 12 

Obviously it is continuous at T c but its derivative is not. 

Now suppose that AM 2 <C /i 2 . Then T c 3> M and so we focus on temperatures 
such that T ^> M. In this case, 

2 / 2/i 2 (l-T/T c ) if T<T C 

™ ef rU J - | ^ (T/Tc _ x) if T > T c l4 '^ 

with T c = Hy/n/j 2 /3AM . The condensate is determined by 8Av 2 (T) = m 2 ff (T) when 
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T <T C while v(T) = when T > T c . The equation of state below T c is 
n ,™s n 2 , f/j, 2 3AM 2 \ „ 3/i 2 MT 



90 V 24 8?r3 / 87r v^ 
tt 2 ^ ^ 2 3AM 2 \ 3// 2 M 

e(T) = 3-T 4 + ^ + — T 2 



and above T c is 



90 V 24 87r3 



vr 2 4 , //i 2 3AM 2 \ 2 /i 4 



P(T) = —T + [—- 
v ; 90 V 24 !6vr 3 J 16A 

^3 ^ 2 3AM 2 \ 
s(T) = 4-T 3 + 2 T 



90 V 24 167r3 , 

tt 2 ^4 ^ 2 3AM 2 \ 2 
e(T) = 3— T + I 777 - — J T + 



90 V 24 167r3 / 16A 

cv(T) = 12^ + 2^- 3 ^V (4.7) 
KV ; 90 V 24 167r3 / 

(Corrections to these formulas are suppressed by the factor exp(— 47r 2 T 2 /M 2 ).) Once 
again, P, s, and e are continuous at T c but cy is not. The most noticeable difference 
is that in the conventional local field theory the effective mass-squared vanishes as 
|T 2 — T 2 |, whereas for the string field theory it vanishes as \T — T c \, although the 
difference is inconsequential in the limit T — > T c . The difference in the exponents is 
seen in the interaction measure too. 



1 M 



_3_P _ I 12 8vr v^FT c 




if T < T c 



^ 1 -— + ' 6^ + 1 if T>T C 

Of course, our calculation is basically a mean-field approximation so the values of 
critical exponents cannot be taken as being very accurate. 

It is instructive to examine the dependence of the discontinuity in the specific 
heat at the critical temperature as a function of T c /M to see the transition from 
conventional local field theory (T c /M < 1) to the SFT limit (T c /M > 1). In those 
two limits the discontinuity can be calculated analytically. 

\ if T c < M 

c v (T c -) - c v (T c +) = XT* { \ / MX 2 (4.9) 



. if T c > M 
8tt 3 V T c ' 
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The discontinuity decreases monotonically with increasing T c /M when measured in 
units of AT C 3 , the only sensible unit for comparison. 

Now we show some full numerical results which do not make any assumption 
about the magnitude of T c /M. Figure 4 shows the dependence of T c /M on the 
variable /i 2 /AM 2 . It changes from the square-root to linear dependence very rapidly 
when T c /M ~ 0.4. Figure 5 shows the dependence of m 2 ff /2/i 2 on T/T c for a value 
of T c /M <C 1 and for a value T c /M ^> 1. The figure clearly shows the quadratice 
dependence for small T c /M versus the linear dependence for large T c /M. Figure 6 
shows the interaction measure, or deviation from the ideal relativistic equation of 
state e = 3P, for both small and large values of T c /M. Finally, Fig. 7 shows the 
discontinuity in the heat capacity at the critical temperature as a function of T c /M; 
it decreases monontically towards zero as T c /M — >■ oo. 




Figure 4: Scaling of the critical temperature with the parameters. The dependence 
changes from square-root to linear around T c /M ~ 0.4. 
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Figure 5: Dependence of the effective mass on temperature for T c /M = 1/100 
(dashed/red) and T c /M = 10 (solid/blue). 

0.3, . . . , 




'5.0 0.5 1.0 1.5 2.0 

T/T c 



Figure 6: The interaction measure as a function of temperature for T c /M = 1/100 
(dashed/red) and T c /M = 10 (solid/blue). 
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0.6 




T c /M 



Figure 7: Discontinuity in the heat capacity at the critical temperature. 



5 Effective Potential 



For a given temperature the field has a stable equilibrium value v(T), as discussed 
and computed in previous sections. If for some reason the average deviates from its 
thermal value by an amount £, there will be a restoring force. This restoring force is 
described by an effective potential £/(£). In essence, this is an expansion away from 
equilibrium states. It is useful in many areas of physics, including statistical physics, 
particle physics, and cosmology. In this section we compute the first few terms in the 
expansion at the 1-loop order. 

Let us define 

= + 0/ (5.1) 

where is a constant field and 0/ is the fluctuation around it whose average value is 
zero. Let us further write 

= v (T)+£. (5.2) 

Here v(T) is the equilibrium value of the condensate at temperature T which was 
previously determined. The £ is a space and time independent (or slowly varying 
on all natural length and time scales) deviation from the equilibrium value that we 
take as a parameter to be varied at will. We will calculate deviations from the 
thermodynamic potential at each temperature as a function of £; this is the effective 
potential (commonly referred to as the effective action in particle physics). The 
Lagrangian can be written as 

C = £ quad + £ int + Cct - V cl (v) - U cl (0 - VK0/) (5.3) 
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where 



-■quad 2 ( 



e _ D/M2 ( D + /i 2) _ 12At; 2_ 5m 2j 0/ 



vv 4 



£ int = -4Au$-A$ 
£ ct = I( 5m 2 + 7 ) 2 

Ki(^) = -\{^ + l)v 2 + \t 

U d (0 = (4At; 2 -/i 2 -7)<+|(l2At; 2 -/i 2 -7)e 2 + 4A< 3 + Ae 4 
V^ f ) = 6A(2< + e 2 )0j + 4A^ 3 / . (5.4) 

The first four terms in the Lagrangian were already introduced and used previously. 
The last two terms vanish when £ = 0. The £/ c i(£) is the contribution to the classical 
potential from the £ field. The V£(0/) gives the interaction between the deviational 
field £ and the quantum field 0/; it will be used to determine the effective potential 
at a given temperature. 

The Feynman rules corresponding to the above action are the same as before but 
with additional terms arising from the presence of £. The thermodynamic potential 
is now 

n = V cl (v) + C/ cI (0 - ^ In y exp Qf " dr d 3 x [£ quad + C mt + C ct - Vfl 

(5.5) 

In a diagrammatic expansion the field 0/ is represented by a solid line while the 
external field £ is represented by a wavy line. The vertices can easily be read off from 
the expressions above. For example, the quartic interaction 0j- has the vertex —A, 
the cubic interaction (j) 3 has the vertex — 4Xv, and the cubic interaction £0 2 has the 
vertex — 12At> . 

Now let us extrapolate away from the equilibrium value of the condensate so that 
£ is not equal to zero. The effective potential for £ is 

U(0 = C/ c i(0 + C/i oop (0 (5-6) 

where 

Uloo » {0 ~ -v ln \ Jrn^s j • ( 5 - 7 ) 

Here S is the action due to £ qua d + Ant + £ct- This can be expanding in an infinite 
series in £. The effective potential has the property that U (0) = 0. 



It is straightforward to compute the 1-loop contribution to U to all orders in £. 
Expand the exponential of to the Aq'th order in £0 2 and to the iVyth order in £ 2 2 ; 
the term £0^ cannot contribute to the one loop order. Expansion of the exponential 
gives rise to a factor l^AqlA^!). Each of these N\ + N 2 terms has two <pf legs. These 
must be connnected to make one and only one loop. The ordering does not matter. 
Taking into account the vertices, this leads to 

{N x + N 2 -1)\. nAX .. Nl/ 10 .^JV 2 ^v^ f d 3 p 



2iVi!iV 2 ! 



24XvO Nl {-12\e) N ' 2 Tj2 [ 4^3 VNl+N2 ( 5 - 8 ) 
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Summing over all N\ and N 2 gives {7i_i oop (0 with the obvious requirement that 
N1 + N2 > 0. When summed with {7 c i(0 the term linear in £ should vanish, otherwise 
we would not be at the extremum of Q. This term is 

which does vanish on account of the equation satisfied by v at one loop order. 

Let us examine U\ oop (^) when T > T c . Since v(T) = there are no three point 
vertices. 

(5.9) 

The reason for separating out the term quadratic in £ is that it naturally combines 
with the quadratice piece in U cl to yield |m 2 ff £ 2 . The remaining terms can be exactly 
summed; they are referred to as the ring diagrams 

d 3 p 



The propagator is 



(2tt)* 



[ln (1 + 12\£ 2 V) - 12\£ 2 V] 



V 



-1 



/ 2 1 2 2\ 1 2,2 



(5.10) 



(5.11) 



The sum and integral in U r i ng are dominated by n = and p — » 0, respectively, 
because there is no need for a UV cut-off in (5.9) for N >2 . The result is 

d 3 p 



U Tmg {^T) = \T J 



(2tt) ; 



In 1 



me 



12Af 



pZ _|_ m z^ j pZ _|_ m 



cff 



12tt 



(m e 2 ff + 12A£ 2 ) 3/2 - m c 3 ff - 18Am eff £ 2 



Thus the potential for T > T c is 

E/(£,T) = ±m 2 ff (T)£ 2 + A£ 4 + *W£,T) + 



(5.12) 



(5.13) 



One must be careful about using this expression too close to T c where fluctuations 
are large and where critical phenomena occur (nonanalytic critical exponents etc.). 
It is only valid when £ 2 < m 2 ff /12A because otherwise the series does not converge to 
a logarithm. This region shrinks to zero as T c is approached from above. 

When T < T c a double series must be summed. In this case we write 

CW£,T) = 6\(2v£ + e)Tj2 J (l^ + fW^T) (5.14) 
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where 



£W«.r) = -ir£ J 



+ E 

JVi=l 



d 3 p 
(2tt) 3 



oo oo 

EE 

,JVi=l 7V 2 =1 



(JVl + Af2 - 1)!( _^ i( _^ 



iVi!iV 2 ! 



+ E 



A ? 2 



1 Af 2 =l 



+ x + y 



(5.15) 



and where x = 24\v£,V and ?/ = 12A£ 2 D. By using the integral representation 

POO 

(Aq + iV 2 - 1)! = / die" 4 ^1+^-! (5.16) 
the sums can be done followed by integation over t with the result that 



(5.17) 



As before, the sum and integral in U T [ ng are dominated by n = and p — » 0, with the 
result that 



T 



12vr L 



(m 2 ff + 24A< + 12A£ 2 ) 



2\3/2 3 



- |m eff (24A< + 12A£ 2 ) 

(5.18) 

Thus the potential for T < T c , including both the classical and one loop contributions, 
is 

U(£, T) = \ml a {T)e + 4A< 3 + A£ 4 + [/ ring (£, T) + • ■ ■ (5.19) 

This expression makes use of the solution for v(T). Previous caution concerning the 
radius of convergence in £ apply here as well. 

Both the T < T c and T > T c potentials have the property that 

d 2 U(i = 0,T) 



which is an oft-cited relationship. However, it is only true when the self-energy is 
frequency and momentum independent. The astute reader will notice that there is 
a contribution of order £ 2 coming from £/ r i ng when T < T c (but not when T > T c .). 
This contribution is of order A 3//2 and therefore is subleading in an expansion in A. 
The magnitude and sign of this contribution is exactly that arising from the one 
loop self-energy diagram with cubic vertices, evaluated in the zero frequency and zero 
momentum limit; see eq. (3.22). 



<9 2 £/ ring (e = 0) 



n 



cubic 



(0,0) 



(5.20) 



We have already decided to drop such subleading terms. For further discussion on 
this point see j38j [36] 
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It should be apparent that the results of this section are independent of whether 
the underlying action is the string field theory or the conventional local field theory. 
The differences only appear when explicit functions of m e s(T) and v(T) are used 
along with the relationship of T c to the paramters in the action. For example, when 
T C <CM (and dropping the ring contribution) 











2 /I 2 


ft 






ft 




1 2 f 





+ V2 



if T < T, 



if T > T r 



(5.21) 



and when T c 3> M (again dropping the ring contribution) 



U 





(1 


2 /i 2 






ft 


2 /i 2 





271-^ T c m eff 
3 M fi 



+ 



7T\/7r T r 



3 M 



3 M 



if T <T C 

if T > T c 
(5.22) 



As examples, we show the effective potential (without the ring contribution) in 
Figs. 8 and 9 for T c /M = 1/100 and T c /M = 10, respectively, for T below, at, and 
above the critical temperature. 



0.15 




Figure 8: Effective potential for T c /M = 1/100 at T/T c =0.5 (solid/blue), 1.0 (dash- 
dotted/black), and 1.5 (dashed/red). 
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Figure 9: Effective potential for T c /M = 10 at T/T c =0.5 (solid/blue), 1.0 (dash- 
dotted/black), and 1.5 (dashed/red). 



6 Conclusion 

We investigated the thermodynamic properties of a nonlocal tachyon in string field 
theory. We first studied the phase transition in the high temperature limit to 1-loop 
order. We introduced the mass shift Sm 2 and calculated the thermodynamic potential 
around the temperature dependent true minimum, v(T). We found that at the 1-loop 
level both the mass shift, 5m 2 , and the minimum, v(T), go to zero at the same critical 
temperature which strongly indicates a second order phase transition. 

We then performed a more involved analysis. Here we included a counter term 
to allow us to make a comparison to the conventional </> 4 theory. We calculated the 
equation of state to 2-loop order. We argued that the cubic contribution is suppressed 
compared to the quartic near the critical temperature. We were able to analytically 
calculate the equation of state and hence the thermodynamic quantities in the limit 
T, M > /i. As expected, we found that at the critical temperature both the pressure 
and the entropy density were continuous but the heat capacity was not, signaling a 
second order phase transition. The discontinuity in the heat capacity, in natural units 
of AT C 3 , was found to decrease monotonically with T c /M . 

Checking the consistency of our analysis, we made the comparison with ordinary 
local field theory by taking the limit that the string parameter M — » oo. Doing so 
we recovered the usual local field theory result. We compared this with analytical 
approximations of our thermodynamic quantities for finite M. We also calculated the 
interaction measure for both cases. Aside from the M dependence of the thermody- 
namic quantities in the string case, we found that the effective mass-squared vanished 
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as \T — T c \ compared to |T 2 — T 2 | in the conventional case. 

In the last section we included the possibility of the the field being out of equilib- 
rium at temperature T by a small amount and computed the corresponding effective 
potential. This allowed us to compute the finite temperature effective potential. We 
found that we were able to calculate the 1-loop contribution to the effective potential 
at all orders in £ for temperature both above and below T c . For both cases the ring 
contribution is only valid for £ 2 < m 2 e ^/\2\. It was also seen that both potentials 
satisfy the usual relationship m 2 ff (T) = d 2 U{^ = 0,T)/<9£ 2 . 

We were able to calculate results that were consistent with conventional scalar 
field theory in the relevant limit. In the limit investigated, T, M > fx, we found that 
this nonlocal theory is very similar to the conventional one, but we were able to see 
effects from the stringy nonlocality. The formalism developed in this paper will help 
us explore the more challenging, but perhaps more interesting, case when M ~ //. 
Our calculations may also be relevant for capturing the thermal properties of the 
Early Universe which is relevant for some cosmological models . 

This work was supported by the U.S. DOE Grant No. DE-FG02-87ER40328. 



References 

A. Sen, Int. J. Mod. Phys. A20, 5513-5656 (2005). 

A. Sen, JHEP 0204, 048 (2002). 
A. Sen, JHEP 0207, 065 (2002). 

E. Witten, Nucl. Phys. B268, 253 (1986). 

V. A. Kostelecky and S. Samuel, Phys. Lett. B207, 169 (1988). 
V. A. Kostelecky and S. Samuel, Nucl. Phys. B336, 263 (1990). 

I. Y. Aref'eva, A. S. Koshelev, D. M. Belov and P. B. Medvedev, Nucl. Phys. 
B638, 3 (2002). 

G. Calcagni and G. Nardelli, Phys. Rev. D 78, 126010 (2008). 

W. Taylor and B. Zwiebach, "D-branes, Tachyons, and String Field Theory," in 
Strings, Branes and Extra Dimensions (TASI 2001), World Scientific, ed. S. S. 
Gubser and J. D. Lykken. 

G. Calcagni and G. Nardelli, Nucl. Phys. B823, 234 (2009). 



B. Dragovich, arXiv:hep-th/0703008. 



B. Dragovich, Theor. Math. Phys. 157, 1671 (2008). 
[10] N. Moeller and B. Zwiebach, JHEP 0210, 034 (2002). 



28 



[11] I. Y. Aref'eva and L.V. Joukovskaya, JHEP 0510, 087 (2005). 

I. Y. Aref'eva, A. S. Koshelev and S. Y. Vernov, Phys. Rev. D 72, 064017 (2005). 
I. Y. Aref'eva, |arXiv:astro-ph/0410443| y2. 

[12] N. Barnaby, T. Biswas and J. M. Cline, JHEP 0704, 056 (2007). 
J. E. Lidsey, Phys. Rev. D 76, 043511(2007). 

N. J. Nunes and D. J. Mulryne, AIP Conf. Proc. 1115, 329 (2009). 
N. Barnaby and J. M. Cline, JCA 0707, 017 (2007). 
N. Barnaby and J. M. Cline, JCA 0806, 030 (2007). 

[13] I. Y. Aref'eva and L. V. Joukovskaya and S. Y. Vernov, JHEP 0707, 087 (2007). 

[14] L. Joukovskaya, Phys. Rev. D 76, 105007 (2007). 
A. Koshelev, JHEP 0704, 029 (2007). 

[15] G. Calcagni, M. Montobbio and G. Nardelli, Phys. Rev. D 76, 126001 (2007). 

[16] G. Calcagni and G. Nardelli, Int. J. Mod. Phys. D 19, 230 (2010). 

[17] I. Y. Aref'eva, L. V. Joukovskaya, and S. Y. Vernov, J. Phys. A: Math. Theor. 
41, 304003 (2008). 

[18] J. W. Moffat, Eur. Phys. J. Plus 126, 53 (2011). 

J. W. Moffat and V. T. Toth, |arXiv:0908.0780| [hep-ph]. 



J. W. Moffat and V. T. Toth, arXiv:0812.1991 [hep-ph]. 



J. W. Moffat, |arXiv : 0709 .4269 [hep-ph]. 
[19] S. Deser and R. P. Woodard, Phys. Rev. Lett. 99, 111301 (2007). 
[20] Tomi Koivisto, Phys. Rev. D 77, 123513 (2008). 



[21] T. Biswas, E. Gerwick, T. Koivisto and A. Mazumdar, arXiv:1110.5249 [gr-qc]. 
T. Biswas, A. Mazumdar and W. Siegel, JCAP 0603, 009 (2006). 

[22] J. Khoury, Phys. Rev. D 76, 123513 (2007). 



[23] E. Elizalde, E. O. Pozdeeva and S. Y. Vernov, [arXiv: 1110 .5806 [astro-ph.CO]. 
[24] J. W. Moffat, Eur. Phys. J. Plus 126, 43 (2011). 



[25] L. Modesto, |a7Xiv: 1107. 2403] [hep-th]. 

[26] S. Capozziello, E. Elizalde, S. 'i. Nojiri and S. D. Odintsov, Phys. Lett. B671, 
193 (2009). 



29 



[27] P. G. O. Freund and M. Olson, Phys. Lett. B199, 186 (1987). 

P. G. O. Freund and E. Witten, Phys. Lett. B199, 191 (1987). 

L. Brekke, P. G. O. Freund, M. Olson and E. Witten, Nucl. Phys. B302, 365 
(1998). 

[28] T. Biswas, J. A. R. Cembranos and J. I. Kapusta JHEP 10, 048 (2010). 

T. Biswas, J. A. R. Cembranos and J. I. Kapusta Phys. Rev. Lett. 104, 021601 
(2010). 

[29] J. J. Atick and E. Witten, Nucl. Phys. B310, 291 (1988). 
[30] J. A. Minahan, JHEP 0103, 028 (2001). 

[31] T. Biswas, M. Grisaru and W. Siegel, Nucl. Phys. B708, 317 (2005). 
[32] Robert Bluhm, Phys. Rev. D 43, 4042 (1991). 

[33] A. Nayeri, R. H. Brandenberger and C. Vafa, Phys. Rev. Lett. 97, 021302 (2006). 

R. H. Brandenberger, S. Kanno, J. Soda, D. A. Easson, J. Khoury, P. Martineau, 
A. Nayeri and S. P. Patil, JCAP 0611, 009 (2006). 

T. Biswas, R. Brandenberger, A. Mazumdar and W. Siegel, JCAP 0712, 011 
(2007). 

J. Magueijo and L. Pogosian, Phys. Rev. D 67, 043518 (2003). 

J. Magueijo and P. Singh, Phys. Rev. D 76, 023510 (2007). 

Y.-F. Cai, W. Xue, R. Brandenberger abd X.-M. Zhang, JCAP 0906, 037 (2009). 

J. Magueijo, L. Smolin and C. R. Contaldi, Class. Quant. Grav. 24, 3691 (2007). 

J.-P. Wu and Y. Ling, Phys. Lett. B684, 177-180 (2010). 

J. Magueijo, Phys. Rev. D 76, 123502 (2007). 

[34] A. Berera, Phys. Rev. Lett. 75, 3218 (1995). 

[35] Y. Volovich, J. Phys. A 36, 8685 (2003). 

V. S. Vladimirov and Y. I. Volovich, Theor. Math. Phys. 138, 297 (2004) [Teor. 
Mat. Fiz. 138, 355 (2004)]. 

V. S. Vladimirov, Izv. Math. 69, 487 (2005). 



D. V. Prokhorenko, arXiv:math-ph/0611068 



N. Barnaby and N. Kamran, JHEP 0812, 022 (2008). 

N. Barnaby and N. Kamran, JHEP 0802, 008 (2008). 

G. Calcagni, M. Montobbio and G. Nardelli, Phys. Lett. B662, 285 (2008). 

P. Gorka, H. Prado and E. G. Reyes, J. Math. Phys. 51, 103512 (2010). 



30 



[36] J. I. Kapusta and C. Gale, Finite Temperature Field Theory, Cambridge Univer- 
sity Press, Cambridge, 2nd edition, 2006. 

[37] T. D. Lee and M. Margulies, Phys. Rev. D 11, 1591 (1975). 

[38] P. Arnold and O. Espinosa, Phys. Rev. D 47, 3546 (1993); Erratum 50, 6662 
(1994). 



31 



